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Abstract 

We present a next-to-leading order determination of the gluino-mediated SUSY 
contributions to Bd - Bd mixing and to the CP asymmetry aj/^x s i R the framework 
of the mass-insertion approximation. Using hadronic matrix elements recently com¬ 
puted on the lattice, we obtain improved constraints on the squark-mass splittings. 


1 Introduction 


With the advent of B factories, B physics is playing a key role in testing the Standard 
Model (SM) picture of flavour and CP violation and in probing virtual effects from 
new physics at low energies. In particular, measurements of A m dl the mass difference 
in the B d - B d system, and of the time-dependent CP asymmetry in Bd —> 

J/i/jK s decays, can be used to put stringent constraints on new-physics contributions 
to A B = 2 processes. In terms of the matrix elements of the effective A B = 2 
Hamiltonian we have 


A m d = 2Abs[(B d |W e A ff B = 2 |B d >] (1) 

a J/^K a = sin 2 f3 eS sin A m d t , (2) 

where 

2p eS = Arg[(B d |7f e AB = 2 |B d )]. (3) 

Within the SM the value of sin 2/3 e ff = sin 2/3 can be connected to the CP violation 
angle <5 of the Unitarity Triangle (UT) 


sin 2/3 


2 a sin <5(1 — a cos <5) 
1 ± a 2 — 2a cos <5 


(4) 


where a = \Jp 2 + rj' 2 in the Wolfenstein parametrisation Q. Assuming that there 
are no new physics effects in B d - B d and K° - K° mixing, sin 2/3 is well determined 
by the standard UT analysis. For example, ref. |2| quotes 


sin 2/3 = 0.698 ±0.066. 


(5) 


The first direct measurements of aj/^K a at the B factories suggested a rather low 
value of the asymmetry || which, if confirmed, would have been a hint of new physics 
0-Q. The most recent world average Q 

sin 2/3 = 0.79 ±0.10, (6) 

however, is in very good agreement with eq. (|5|) and therefore does not favour extra 
contributions. In spite of this, it is worth investigating, in neutral B meson systems, 
signals of new physics which may emerge in the near future thanks to decreasing 
experimental errors and theoretical uncertainties. On the one hand, for A B = 2 
transitions, where the SM seems to agree with the data, it is important to improve 
the theoretical accuracy for constraining new physics contributions to box diagrams. 
On the other, more accurate predictions of mixing help to uncover signals of new 
physics in B decay asymmetries, for example B —> Kit or B —» K<p. These processes 
are dominated by loop contributions and one measures simultaneously the effect 
of mixing and CP violation in decay amplitudes, from both SM and new physics 
virtual particles. Finally, of particular interest is the B s -B s mixing for which the 
SM amplitude is real and the phase entirely originates from new physics effects. This 
may give rise to CP violating effects which would be absent in the SM. 

In this work, we aim at testing new-physics effects in A B = 2 transitions using the 
most recent experimental and theoretical results: the new data from B factories and 
the very recent lattice determinations of the relevant matrix elements, which allow 
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the computation of Am^ and of aj/^K a at the Next-to-Leading Order (NLO), apart 
from some additional model-dependent 0(a s (M\y)) terms. First of all, we present 
model-independent formulae for A B = 2 ( B c i and B s ) transitions which include all 
possible new-physics contributions. These formulae are obtained by considering the 
most general effective Hamiltonian for A B = 2 processes, 7i^ B=2 . We write down 
the expression for Am^ (A m s ) and for /3 e g in terms of the Wilson coefficients at 
the electroweak scale, including NLO QCD corrections in the running from M-\y to 
H = TTZfo ]|, |^] and the matrix elements very recently computed in lattice QCD [ 101, 
in the same RI-MOM renormalisation scheme. These expressions can be readily used 
to compute Aand /3 e g in any extension of the SM. 

As a concrete example of physics beyond the SM, we then consider the Minimal 
Supersymmetric Standard Model (MSSM) with arbitrary soft breaking terms and 
we study model-independent constraints on the SUSY parameter space coming from 
A md and aj/^x a - We closely follow the analysis of ref. [11], where upper bounds 
on new sources of flavour and CP violation were obtained from the study of K°— 
K° mixing. For earlier analyses, either at tree level or with LO evolution, using 
naive Vacuum Insertion Approximation (VIA) B parameters, see refs. 112] and [13] 
respectively. Ref. [11] and the present work can be considered as first steps toward 
a full NLO model-independent analysis of flavour and CP violation beyond the SM, 
with special focus on its most studied extension, namely the general MSSM. We 
stress that the inclusion of hadronic matrix elements from lattice QCD, for operators 
renormalised consistently with the Wilson coefficients at the NLO, is highly important 
to reduce the uncertainties in the study of flavour and CP violation. Indeed, while 
matrix elements for operators appearing in the SM have been studied for a long time, 
it was only recently that similar results have been obtained for operators that only 
come from physics beyond the SM [1C, |h|]. More work is still to be done along these 
lines. For example, the large uncertainties in the SUSY constraints from the electric 
dipole moment of the neutron could be drastically reduced by a lattice computation 
of the relevant matrix elements [ 151. 

The paper is organised as follows. In the next section we introduce the most 
general 7i^ =2 with the inclusion of the NLO QCD corrections in the evolution 
from the scale of new physics down to low energy. We provide the expression of 
the Wilson coefficients at the low hadronic scale (~ mb) as a function of the Wilson 
coefficients and of a s at the “new physics scale”. Section || deals with the evalu¬ 
ation of the hadronic matrix elements of the local operators appearing in 7i^ =2 . 
Here we replace the “traditional” values of the B parameters in the VIA with the 
values recently obtained from a lattice computation in the RI-MOM scheme [lOj. 
We then turn to SUSY and present in Section f| the expression of gluino-mediated 
contributions to 7i^ =2 . As long as we deal with general squark mass matrices, the 
inclusion of the gluino-mediated FCNC alone is sufficient to get the correct bulk of 
the SUSY contribution. In specific models corresponding to particularly restricted 
squark mass matrices, however, there are regions of the SUSY parameter space where 
charginos, stops and/or charged Higgs are relatively light and their exchange gives 
contributions even larger than gluinos to FCNC |n|. This happens more easily in 
the case of B^ - B^ mixing, rather than for K°-K°, because the relevant mixing 
angles and/or Yukawa couplings are larger in this case. Unfortunately, the experi¬ 
mental information and the theoretical accuracy are presently insufficient to attempt 
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a simultaneous analysis which includes a spanning of the squark mass matrices and 
of the stop/chargino mass values. For this reason, in this first study which uses the 
new measurements of sin 2/3, NLO corrections and lattice B parameters, we include 
gluino-mediated FCNC only. Our quantitative results are presented in Sect. |] and 
summarised in Tables ^ and [| We perform our computations in the mass insertion 
approximation [[Id]]. We choose the super-CKM basis for the fermion and sfermion 
states, where all the couplings of these particles to neutral gauginos are flavour di¬ 
agonal, while genuine SUSY FC effects are exhibited by the non-diagonal entries of 
the sfermion mass matrices. Denoting by A 2 the off-diagonal terms in the sfermion 
mass matrices (i.e. the mass terms relating sfermions of the same electric charge, 
but different flavour), the sfermion propagators can be expanded as a series in terms 
of 5 = A 2 /m 2 , where rh is the average sfermion mass. As long as A 2 is signifi¬ 
cantly smaller than m 2 , we can just take the first term of this expansion and then 
the experimental information concerning FCNC and CP violating phenomena trans¬ 
lates into upper bounds on the 5s [11, 12]. Tables [2] and § contain new constraints 
on the parameters 5. With respect to previous analyses, we find that the inclusion 
of the direct measurements of sin 2 (3 allows imposing quite stringent constraints on 
Im5i3 [ 181 - On the theoretical side, as already observed in [111, the QCD evolution 
from large scales and the use of the lattice B parameters induces sizable changes in 
the limits on the values of the 5s. The major effect of our improved computation is 
felt by Left-Right operators. 


2 Effective Hamiltonian for A B — 2 processes 
beyond the SM 

The most general effective Hamiltonian for A B = 2 processes beyond the SM has 
the following form: 

5 3 

H * S B = 2 = X) C * Qi + Ci Qi (?) 

i=l i =1 

where 


Qi 

Q'2 

Qs 

Qa 

Q5 


dl^btdir^L , 


4444 , 
4444, 
4444> 
4444, 


( 8 ) 


and the operators Q 1,2,3 are obtained from the Qi, 2,3 by the exchange L <-> R. Here 
Qr,l = Pr.l q, with Pn t L = (1 ± 75)/2, and a and f3 are colour indexes. In the case 
of the B s system one has simply to replace the d with the s quark in the operators 
appearing in eqs. (|8|). 

The NLO anomalous dimension matrix for the most general Tffg' =2 has been 
computed in J8| and the results have been confirmed in ||. We use the Regularisation- 
Independent anomalous dimension matrix in the Landau gauge (also known as RI- 
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Constant 

Value 

Constant 

Value 

m Bd 

5.279 GeV 

fs d 

200 ± 30 MeV 

mfS( 2GeV) 

7 MeV 

m MS( m MS) 

4.23 GeV 

m t 

174 ± 5 GeV 

Oi s {.M z ) 

0.119 

ivy 

(40.7 ± 1.9) x 10" 3 

\Vub\ 

(3.61 ±0.46) x 10~ 3 


Tabic 1: Input quantities used in the phenomenological analysis. m^ IS and m^ IS are the 
MS masses whereas rrit is the pole mass. 


MOM), since we will make use of matrix elements computed in lattice QCD with the 
same choice of renormalisation scheme. 

A full NLO computation would also require the 0(a s ) corrections to the matching 
conditions which determine the Wilson coefficients, see eqs. (15) below. These correc¬ 
tions are model-dependent and have been computed only in few specific cases |l9|, ^Of . 
One might argue that, being of order a s {Mg) (here and in the following Ms repre¬ 
sents the scale of new physics), these contributions should be small, as suggested by 


the cases of the SM and of the two Higgs doublet model [19]. This statement, how¬ 
ever, can only be confirmed by explicit computations in specific models. Due to the 
absence of 0(a s ) corrections to the matching, our Ti^ =2 are affected by a residual 
scheme dependence, which would be cancelled by the terms of order a s {Ms ) in the 
Wilson coefficients Ci(Ms). These terms are model dependent. 

The Ci(Ms) are obtained by integrating out all new particles simultaneously 
at the scale Ms- We then have to evolve the coefficients down to the hadronic 
scale p = mb = 4.6 GeV (mb = mb(p = mb) is the RI-MOM mass), which is 
the renormalisation scale of the operators used in ref. [1C ]. By coincidence, mb 
has the same numerical value as the pole mass, m^ olc = 4.6 GeV, extracted at the 
NLO from the MS mass m^ IS (m^ IS ) = 4.23 GeV, see Table 0 For consistency, we 
have to evolve the Wilson coefficients at the same renormalisation scale. The SM 
contribution can be computed independently and evolved from M\y to p using the 
well-known NLO QCD corrections |H]]. 

We give here an analytic formula for the expression of the Wilson coefficients at 
the scale p = mb as a function of the initial conditions at the SUSY scale Ci(Ms) 
and of a s (Ms ). This formula has been obtained by using, for the SM parameters, 
the values in Table |I[ For Ms > mt we obtain 


C r (m p b ole ) = ( 6 S r ’ S) + hcp s) ) l ai C S (M S ), 


(9) 


where, in the evolution of the coefficients from Ms, we have chosen Ms = (Mg+Mq )/2 
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and r] = a s (Mg)/a s (mt)■ The magic numbers are given below: 
at = (0.286,-0.692,0.787,-1.143,0.143) 


b\ U) = (0.865,0,0,0,0), 

c( n) = (-0.017,0,0,0,0), 

bf 2) = (0,1.879,0.012,0,0), 

cf 2) = (0,-0.18,-0.003,0,0), 

&f 3) = (0,-0.493,0.18,0,0), 

cf 3) = (0,-0.014,0.008,0,0), 

bf 2) = (0,-0.044,0.035,0,0), 

cf 2) = (0,0.005,-0.012,0,0), 

bf z) = (0,0.011,0.54,0,0), 

cf 3) = (0,0.000,0.028,0,0), 

^ 44) = (0,0,0,2.87,0), 

c ( 44) = (0,0,0,-0.48,0.005), 

b[ i5) = (0,0,0,0.961,-0.22), 

c ( 45) = (0, 0,0, -0.25, -0.006), 

bf 4) = (0,0,0,0.09,0), 

cf 4) = (0, 0,0, -0.013, -0.016) 

b[ 55) = (0,0,0,0.029,0.863), 

cf 5) = (0,0,0,-0.007,0.019), 


( 10 ) 


and we have only written the non-vanishing entries. The magic numbers for the 
evolution of C 1-3 are the same as the ones for the evolution of Ci_ 3 . Formulae (j£) 
and (1C) are combined with the B parameters evaluated at m^ ole = 4.6 GeV in 
the MOM-RI scheme [1C], given in eqs. (jb|) below 1 . In this way it is possible to 
compute the contribution to A md, s and /3 e ff at the NLO in QCD for any model of 
new physics in which the new contributions with respect to the SM originate from 
extra heavy particles. One just has to plug in the expression for the C) evaluated 
at the large energy scale Mg in his favourite model. When the 0(a s ) corrections 
to the Ci(Mg) are available, one can obtain a full NLO, regularisation-independent 
result; in the cases where this corrections have not been computed yet, the results 
contain a residual systematic uncertainty of order a s (Adg). We note that, due to 
the presence of large entries in the NLO anomalous dimension matrix, a systematic 
uncertainty is present in the QCD evolution from the SUSY scale to the hadronic 
one. This uncertainty stems from different ways of writing the NLO evolution matrix, 
which differ by terms of O(o^), and has been taken into account in the numerical 
analysis as uncertainty in the numerical value of the magic numbers. For this reason 
we have presented only the significant figures of the magic numbers, which have an 
uncertainty of one unit on the last digit. The only exception is c^ 11 ^ for which the 
uncertainty is 3 units. 


3 Hadronic Matrix Elements 

The matrix elements of the operators Qi between neutral B mesons in the VIA are 
given by: 

1 


{Bd\Qi\B d )yiA = 7; m B d fB d , 


(Bd\Q2\Bd)viA — 


™B d 


24 \ nit + m d 


{Bd\Q3\Bd)viA = 


1 




_ 24 \mb + m d 

1 Magic numbers in the NDR scheme can be found in ref. 


m B d f 2 B d , 


m B d fl d , 


( 11 ) 
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{Bd\Q4\Bd)viA 

{Bd\Q$\Bd)\riA 


J_ + I f m B d 

24 4 \rrib + md 


2 


-+—( mBd y 
8 12 / 


m B d f 2 Bd , 


m B d fl d , 


where is the mass of the meson and m b and md are the masses of the b 
and d quarks respectively. Note that the normalisation of the operators used in this 
paper differs by a factor of 2 ms d from the one used in ref. 11C ]. Here and in the 
following, the same expressions of the B parameters of the operators Q 1-3 are valid 
for the operators Qi~ 3 , since strong interactions preserve parity. In the case of the 
B s system one has simply to replace the d with the s quark in the operators and 
expressions appearing in eqs. ([TT|) . 

In the case of the renornralised operators, we define the B parameters as follows: 


(BdlQ^lBd) 
(B d \Q 2 (li)\B d ) 
(.B d \Q 3 (ii)\B d ) 
{B d \Q^)\B d ) 
(B d \Q 5 (ri\B d ) 


m B d fB d Bl{n) , 


_5_ 

"24 




1 

24 

1 




4 \m b (n) + m d {n) 


m B d fB d B 4 ( m ) 


12 




( 12 ) 


where the notation Qi(/r) (or simply Qi) denotes the operators renormalised at the 
scale /r. For consistency with the calculation of the B parameters in ref. [1C], the 
quark masses, evolved at the scale //, must be computed in the Landau RI-MOM 
scheme p3[ p4| . With the numbers in Table II this corresponds to m b = mb(m b ) = 
4.6 GeV and m d (m b ) = 5.4 MeV (RI-MOM) 0 

A few words of explanation are necessary at this point. The B parameter of 
the matrix element {B d \Qi\Bd) has been extensively studied on the lattice due to 
its phenomenological relevance 25, 26], and used in many UT analyses [0, 27]. For 


the other operators, instead, all the studies beyond the SM have used the VIA B 
parameters, which in some cases, as will be shown below, is a crude approximation. 
The numerical results for the B parameters Bj(^) refer to the RI-MOM scheme. We 
have used the same definition of the B parameters as in ref. [10] from which we have 
taken the values. In our numerical study, we have used 


Bi(m b ) = 0.87(4), B 2 {m b ) = 0.82(3)(4) , 

B 3 {m b ) = 1.02(6)(9) , B 4 {m b ) = 1.16(3)1® , (13) 

B 5 {m b ) = 1.91(4)lf . 


We added a conservative ±10% systematic error to the errors quoted above, based 
on previous experience on discretisation, extrapolation and quenching effects. In 

2 The precise value of is not very important for our analysis. 
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the cases where a comparison has been possible, the quenching effects for all the 
B parameters computed so far have always been found much smaller than all other 
uncertainties. 


4 Effective Hamiltonian for A B = 2 processes 
in SUSY 


In this Section, we describe the computation of gluino-mediated contributions to 
•?A/= 2 at the NLO in QCD. 

Several cases may be considered: i) rrig ~ rrig, ii) rrig <C rrig and iii) rriq 3> nig , 
where rriq is the average squark mass and nig is the gluino mass. Case ii), in which 
rrig <C nig, cannot be realised in the framework considered here, in which the soft 
SUSY breaking terms are introduced at the Planck scale. This is due to the fact that 
the evolution from the Planck to the electroweak scale forbids such a mass hierarchy. 
In fact, neglecting the effects of Yukawa couplings and A-terms, one obtains in the 
down sector the following approximate expression for the ratio x = m~/m~ at the 


electroweak scale, in terms of the value xq at the Planck scale [28]: 


x ~ 


9x 0 


1 - f - 7xq 


9 

7 


(14) 


Even if one starts at the super-large scale with an extreme hierarchy between squark 
and gluino masses (xo 1), at the electroweak scale the two masses will be of the 
same order. For this reason, we will not consider the case rriq rrig in our analysis. 
Case iii), in which rriq S> nig, can be realised in some special models, such as effective 
supersymmetry |3(| or models with a light gluino [31]. The NLO corrections in this 
case have been computed in ref. In this first study, however, we do not consider 
these cases and present only results for the case rriq ~ rrig. 

The mass insertion approximation, which we have adopted in our analysis, presents 
the major advantage that one does not need the full diagonalisation of the sfermion 
mass matrices to perform a test of FCNC in the SUSY model under consideration. 
It is enough to compute ratios of the off-diagonal over the diagonal entries of the 
sfermion mass matrices and compare the results with the general bounds on the 5s 
which we provide here from A mB d and ajj^K s - This formulation of the mass inser¬ 
tion approximation in terms of the parameters 5 is particularly suitable for model- 
independent analyses, but involves two further assumptions: the smallness of the 
off-diagonal mass terms with respect to the diagonal ones, and the degeneracy of the 
diagonal mass terms in the super-CKM basis. The latter assumption, related to the 
use of the average squark mass fh, is well justified in our case, since, on quite general 
grounds, one does not expect a sizable non-degeneracy of the down-type squarks. In 
SUSY GUT it may happen that the sbottom mass is smaller than the d and s squark 
masses. Also in this case, namely for non-degenerate diagonal mass terms, however, 


it is possible to define a generalised mass insertion approximation [ 291 and apply the 
methods exposed in this paper. 

There exist four different 6 mass-insertions connecting flavours d and b along a 
sfermion propagator: LL , (^13 ) rri (^13 ) lr anc ^ rl’ ^ ie ibexes L an d R 

refer to the helicity of the fermion partners. The amplitude for A B = 2 transitions 
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a) 





c) 



b) 


d) 


Figure 1: Feynman diagrams for A B = 2 transitions, with h, k, l, m 


in the full theory at the LO is given by the computation of the diagrams 
At the lowest order in QCD, in the basis of eq. (|8|), the Wilson coefficients 

by 0: 


Cl = 

216m? 

^24 x/ 6 (x) 

+ 66 / 6 ( 

i) ) (+)L ■ 


C 2 = 

- ^-9 204 x (x) 

216m? ' J 

(«13 

) 2 
) RL 

5 


C 3 = 

36 

216m? 

x h{x) 

,\ 2 
13 ^l ’ 



c 4 = 

a 2 s 

216 m? 

(504x fe( 

x) — 

72/ 6 (x)) (^ 3 ) ll 



-132/e (a: 

) (sis) 


) 





'rlJ 



C 5 = 

a 2 

216m? 

(24 x/ 6 (x 

) + 120 fe(x)) (sfs) ll 

(^ 3 ) RR 


—180 fe(x 

>(+L 

« 

) iJL 

5 



61 = “21S^( 24l/ « (l)+66/6W )(' 5 f3)^, 
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in fig. [1|. 
are given 







































0< y< 90 
90 < y < 180 
180 < Y^ 270 
270 < Y< 360 


0.01 


Re ^ 13^LL=RR 



Figure 2: Allowed regions in the ( 7 , Re(Sf :i ) sl, Im(5f 3 ) ll) space with ( 5 f 3 )i,L only (left) 
and (^ 13 )ll = (Sf 3 )rr (right). The two lower plots are the corresponding projections in 
the Re(5f 3 )LL~I m (fi 13 ) ll plane. Different colours denote values of 7 belonging to different 
quadrants. 


C 2 

C 3 


~2ffnl 2MXMX)(l3 

- ^,36x/«(x)(«f 3 


LR 


LR 


where x 


m|/m? and the functions fs(x) and fe(x) are given by: 


Mx) 
h (x) 


6(1 + 3x) In x + x 3 — 9x 2 — 9x + 17 
6(x — l) 5 

6x(l + x) In x — x 3 — 9x 2 + 9x + 1 
3(x - l) 5 ' 


(15) 


(16) 


In the absence of 0(ct s ) corrections to the matching, we interpret the C* given above 
as coefficients computed at the large energy scale Ms ~ rriq ~ rrig, i.e. C* = Ci(M$ ). 


5 Numerical Analysis of A B = 2 Processes 

In this section, we illustrate the procedure followed in our analysis and give the main 
results. Arrid and sin 2/3 are defined in terms of the matrix element {B c f\'H / ^ = 2 \B ( i) 
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Figure 3: Allowed regions in the ( 7 , Re(Si 3 ) B R, Im(Sf 3 ) lr) space with (5( 3 )lr only (left) 
and (5f 3 )LR = (Sf 3 )rl (right). The two lower plots are the corresponding projections in 
the Re(5f 3 )LR-Im(5f 3 ) lr plane. Different colours denote values of 7 belonging to different 
quadrants. 

which we schematically write as 

{Rd\RtfP 2 1 Bd) = Re Asm + * Im Asm + 

Asusy Re(df 3 ) 2 AB + i Asusy Im ( 8 f 3 ) 2 AB , (17) 

where Asm is the Standard Model contribution which depends on the CKM matrix 
parameters, in particular on the CP violation phase 7 . (6f 3 ) AB denotes the generic ef¬ 
fective coupling appearing in eq. ©• In the following, for simplicity, the constraints 
are obtained by imposing that the sum of the SUSY contributions proportional to 
a single (Sf 3 ) A B parameter and the SM contribution does not exceed by more than 
la the experimental value of A m,d and sin 2/3. This is justified (a posteriori) by 
noting that the constraints on different (5i 3 )ab parameters exhibit a hierarchical 
structure, and therefore interference effects between different contributions would re¬ 
quire a large amount of fine tuning. We consider the following four cases: i) (Sf 3 )LL 
only; ii) (6( 3 )ll = (Sf 3 )rr\ iii) (Sf 3 ) LR only; iv) (Sf 3 ) LR = ( 6 i 3 )rl- For this reason, 
in all the cases considered here, the amplitude in eq. ( 0 ) only depends on a single 
(^13) AS ’ 

The procedure used in our numerical analysis is the following: 

• We scan over all possible values of 7 between 0 and 2vr since, in the presence 
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2.0 x 10- 2 

1.3 x 10~ 2 

1.2 x 10- 2 


Re(7 3 ) 1H | 

| 

Re(S d i:i ) lr=rl\ 

X 

TREE 

LO 

NLO 

TREE 

LO 

NLO 

0.25 

3.4 x 10~ 2 

2.7 x 10- 2 

3.0 x 10" 2 

3.8 x 10” 2 

2.7 x 10” 2 

2.6 x 10- 2 

1.0 

3.9 x 10" 2 

3.0 x 10“ 2 

3.3 x 10“ 2 

8.3 x 10~ 2 

5.4 x 10~ 2 

5.2 x 10” 2 

4.0 

5.3 x 10“ 2 

4.1 x IQ’ 2 

4.5 x 10“ 2 

1.2 x 10- 1 

2.5 x 10" 1 

— 


1 Im(^ 3 ) LR | 

| 

Im(S 33 ) lr=rl\ 

X 

TREE 

LO 

NLO 

TREE 

LO 

NLO 

0.25 

7.6 x 10" 2 

6.0 x 10~ 2 

6.6 x 10" 2 

1.5 x 10~ 2 

9.0 x 10“ 3 

9.0 x 10~ 3 

1.0 

8.7 x 10" 2 

6.6 x 10- 2 

7.4 x 10- 2 

3.6 x 10- 2 

2.4 x 10“ 2 

2.3 x 10- 2 

4.0 

1.2 x IQ” 1 

9.2 x 10" 2 

1.0 x 10” 1 

2.7 x 10” 1 

5.7 x 10" 1 

— 


Table 2: Maximum allowed values for \ Re (fitj) AB I and \ Ini (8fj) AB |, with A, B = (L, R ), 
for an average squark mass rriq = 500 GeV and for different values of x = We give 

the results in the following cases: i) with the tree level Wilson coefficients, namely without 
evolution from Ms to mb, and VIA B parameters, denoted by TREE; ii) with LO evolution 
and VIA B parameters, denoted by LO; Hi) with NLO evolution and lattice B parameters, 
denoted by NLO. The missing entries correspond to cases in which no constraint was found 

for I {K)ab I < °- 9 ' 

of SUSY effects, this parameter is no longer constrained by the UT analysis. 
Indeed, ex may be subject to large SUSY contributions due to the (Sf 2 )AB 
couplings and the value of 7 can be very different from the SM value; 

• for a fixed value of 7 , we scan over Re(df 3 ) A B and Im(<5 33 )AB between —1 and 1 ; 

• we let the other experimental and theoretical parameters, in particular the B 
parameters, vary within the ranges given in Table |l] and in Eq. (|l3|); 

• in the ( 7 , Re(df 3 ) A b, Im(<ff 3 )AB) space, we select the region where the predicted 
values of Aand sin 2/3 lie within 1 a from the measured values. 

As can be seen from Eq. (0), the amplitude only depends on the square of the 
coupling, (Sf 3 ) AB . Therefore, the selected values of (Sf 3 ) A B have a twofold ambiguity 
(Re(5f 3 ) A B —* > — R e(5i 3 ) A B and Im(d 33 )^B —► — Im((5 33 )^4B) as will be evident from 
our results. 

We are ready to discuss the physics results. In Figs. [ 2 ] ||, for rriq = nig = 500 GeV, 


11 





































1 HeWaM 


X 

rriq = 250 GeV 

rriq = 1000 GeV 

rriq = 250 GeV 

rriq = 1000 GeV 

0.25 

— 

1.3 x 10’ 1 

— 

3.8 x 10~ 2 

1.0 

6.5 x 10" 2 

3.1 x 10- 1 

1.0 x 10~ 2 

4.2 x 10~ 2 

4.0 

3.2 x 10" 1 

— 

1.4 x 10~ 2 

5.9 x 10~ 2 


1 —1 

5^ 

0 

tV 

I m (<5f 3 )LL=RR 

X 

nriq = 250 GeV 

rriq = 1000 GeV 

rriq = 250 GeV 

rriq = 1000 GeV 

0.25 

— 

2.8 x 10- 1 

— 

1.6 x 10~ 2 

1.0 

1.3 x 10" 1 

5.0 x 10” 1 

4.0 x 10" 3 

1.8 x 10~ 2 

4.0 

1.5 x 10" 1 

— 

6.0 x 10“ 3 

2.5 x 10" 2 


[RepSid 

R e (<^ i 3 )lR=Rl\ 

X 

rriq = 250 GeV 

m,q = 1000 GeV 

rriq = 250 GeV 

rriq = 1000 GeV 

0.25 

— 

6.2 x 10- 2 

— 

5.3 x 10~ 2 

1.0 

1.5 x 10~ 2 

6.8 x 10“ 2 

2.5 x 10“ 2 

1.1 x 10" 1 

4.0 

2.2 x 10" 2 

9.4 x 10“ 2 

— 

— 


1 Im(5f 3 ) LR | 

I —1 

0 
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X 

rriq = 250 GeV 

rriq = 1000 GeV 

rriq = 250 GeV 

rriq = 1000 GeV 

0.25 

— 

1.4 x 10” 1 

— 

2.1 x 10~ 2 

1.0 

3.3 x 10~ 2 

1.5 x 10- 1 

9.0 x 10“ 3 

4.5 x 10~ 2 

4.0 

4.8 x 10" 2 

2.2 x 10' 1 

— 

— 


Table 3: Maximum allowed values for \ Re (Sfj) AR \ and \ Im (8fj) AB \, with A, B = ( L , R ), 
for an average squark mass rriq = 250 GeV or 1000 GeV and for different values of x = 
m?/ml. We give the results with NLO evolution and lattice B parameters. The missing 
entries correspond to cases in which no constraint was found for | (Sff) AR \ < 0.9. 


we show the allowed regions in the ( 7 , Re(df 3 )ABi Im((5f 3 )AB) space and the corre¬ 
sponding projections on the (Re(5f 3 )AB, Im^g)^#) plane. For the readers’ con¬ 
venience, we have coloured regions corresponding to 7 in the four quadrants with 
different colours. An alternative way of presenting our results is given in Fig. |I| 
where the allowed values of Ahs(6 33 )AB are plotted as a function of Arg(df 3 )AB ■ 

We give the maximum allowed values of | Re(<5f 3 )AB| and | I]ti(<5i 3 )ab| in Table [ 2 ] 
for x = 0.25, 1 and 4. We also give results for rriq = 250 GeV and 1000 GeV in Table[3]. 
We find that the constraints on the (Sf 3 ) a b parameters are less effective by a factor 
from 10 to 100 with respect to those extracted on (6f 2 )AB from AS = 2 transitions. 
The reason for this difference is twofold. On the one hand, one naturally expects 
constraints of the order of sin 3 9 C for Bj - B^ mixing and sin 5 9 C for K°-K° mixing, 
dictated by the top contribution in the SM amplitude. On the other, in the K°-K° 
case, LR contributions are chirally enhanced, so that limits coming from chirality¬ 
changing operators are more effective. Obviously, specific models of flavour, where 
(Sf 2 )ab and (<5i 3 )ab are correlated, can constrain the squark-mass splittings more 
severely. In these cases one could even attempt the inclusion of chargino contributions 
to the mixing [18, 33]. Our general formulae for the effective Hamiltonian and the 
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magic numbers remain valid, and the phenomenological analysis can be carried out 
along the same lines. 

We now discuss the difference between our results and those obtained with tree 
level or LO Wilson coefficients and VIA matrix elements. In this case, with our 
definition of the B parameters given in eq. (|T^), with k = ( rrib/mB d ) 2 = 0.76, we 
have used 

B\ /ia = 1.0 , B}{ ta = 1.0 , 

B\ ia = 1.0 , B\ ia = 1.13 , (18) 

Bl IA = 2.14 . 

With NLO corrections and lattice B parameters, the constraints are generically 
looser than using LO and VIA B parameters. The theoretical uncertainties of the 
NLO constraints are, however, much smaller. In the LL=RR case, the NLO bounds 
are even more stringent than the LO ones, in spite of the uncertainties on the B 
parameters. The most striking effect is obtained in the LR=RL case for x = 4. Here, 
the strong cancellation taking place at the LO is exacerbated by NLO corrections, 
jeopardising any possibility to constrain {Sf 3 )LR in this case. Clearly, given the 
cancellation, this result should be interpreted with care and verified by calculations 
in any given model. 

With the future measurements of A m s at proton colliders, and hopefully of the 
corresponding CP violating effects in B s mixing and decays, it will be soon possible to 
perform a similar analysis for the B s system. Work along this line is in progress |h| . 


6 Conclusions 


In this work we have provided an improved determination of the gluino-mediated 
SUSY contributions to Bd - Bd mixing and to the CP asymmetry aj/^K g in the 
framework of the mass-insertion method. The improvement consists in introducing 
the NLO QCD corrections to 7i^ =2 jsl p] and in replacing the VIA matrix elements 


with their recent lattice computation [10]. As a glimpse at Table |] readily reveals, 
these improvements affect previous results in a different way, according to the op¬ 
erators of 7i that one considers. The effect is particularly large for Left-Right 
operators. 


We have provided an analytic formula for the most general low-energy TL^ =2 at 
the NLO, in terms of the Wilson coefficients at the high energy scale. This formula 
can be readily used to compute Aand sin 2/3 in any extension of the SM with new 
heavy particles. 

FCNC and CP violating phenomena (in particular in B physics) are promising 
candidates for some indirect SUSY signal before LHC, and are in many ways com¬ 
plementary to direct SUSY searches. From this point of view the theoretical effort 
to improve as much as possible our precision on FCNC computations in a SUSY 
model-independent framework is certainly worth and, hopefully, rewarding. 
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A bs(8 13 ) LL=RR Abs(5 13 ) LL 



Figure 4: Allowed values of Abs(5f 3 )LR as a function of Arg(5f 3 )AB in the four considered 
cases. Different colours denote values of 7 belonging to different quadrants. 
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